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ABSTRACT 

We investigate the extinction curves of young galaxies in which dust is supphed from Type 
II supernovae (SNe II) and/or pair instability supernovae (PISNe). We adopt Nozawa et al. 
(2003) for compositions and size distribution of grains formed in SNe II and PISNe. We find 
that the extinction curve is quite sensitive to internal metal mixing in supernovae (SNe). The 
extinction curves predicted from the mixed SNe are dominated by Si02 and is characterised 
by steep rise from infrared to ultraviolet (UV). The dust from unmixed SNe shows shallower 
extinction curve, because of the contribution from large-sized 0.1 /im) Si grains. However, 
the progenitor mass is important in unmixed SNe II: If the progenitor mass is smaller than 
~ 20 Mq, the extinction curve is flat in UV; otherwise, the extinction curve rises toward the 
short wavelength. The extinction curve observed in a high-redshift quasar {z = 6.2) favours 
the dust production by unmixed SNe II. We also provide some useful observational quantities, 
so that our model might be compared with future high-z extinction curves. 

Key words: dust, extinction — galaxies: evolution — galaxies: high-redshift — galaxies: 
ISM — supernovae: general — quasars: individual: SDSS J104845.05H-463718.3 



1 INTRODUCTION 

Dust grains play a crucial role in the formation and evolution 
of galaxies. Dust grains absorb stellar light and reemit it in far 
infrared (FIR), controlling the energy balance in the interstellar 
medium (ISM) and protostellar gas clouds. Also, the surface of 
dust grains is a site for an effic ient formation of H2 molecules (e.g. 
[Cazaux & Tiel ensl2002lEooi . which act as an effective coolant in 
metal-poor ISM. Those effects of dust turn on even in a very metal- 
poor environment (~ 1% of the solar metallicity; Hirashita & Fer- 
rara 2002; Morgan & Edmunds 2003), and it is argued that the star 
formation rate is enhanced because of the first dust enrichment in 
the history of galaxy evolution. The first source of dust in the uni- 
verse is Type II supernovae (SNe II) or pair instability supernovae 
(PISNe), since the lifetime of their progenitors is short (~ 10'^ yr). 
In the local universe, dust grains are also produced by evolved low 
mass stars (Gehrz 1989), but this production mechanism requires 
much longer (> 1 Gyr) timescales. The first dust supplied by SNe 



* E-mail: hirashita@u. phys.nagoya-u.ac.jp 

t Postdoctoral Fellow of the Japan Society for the Promotion of Science 
(JSPS). 

X Postdoctoral Fellow of the JSPS for Research Abroad. 



II or PISNe may trigger the formation of low-mass stars (Schneider 
et al. 2003). 

Since dust is important even in the early stage of galaxy evolu- 
tion, it is crucial to know how much dust and what species of dust 
form. Todini & Ferrara (2001), following the method by Kozasa, 
Hasegawa, & Nomoto (1989, 1991), show that dust mass produced 
by a SN II is roughly 0.1-0.4 Mq. They also find that SNe II 
form amorphous carbon with size around 300 A and silicate grains 
around 10-20 A. Schneider, Ferrara, & Salvaterra (2004) extend the 
progenitor mass range to the regime of PISNe (140-260 Mq) and 
find that 30-60 Mq of dust forms per PISN. The grain radii are dis- 
tributed from 0.001 to 0.3 fj,m, depending on the species. The moti- 
vation for considering PISNe comes from some evidence indicating 
that the stars formed from metal-free gas. Population III (PopIII) 
stars, are very massive with a characteristic mass of a few hundred 
solar masses (Bromm & Larson 2004, and references therein). Such 
massive stars are considered to begin pair creation of electron and 
positron after the helium burning phase, and finally end their lives 
with explosive nuclear reaction disrupting the whole stars (Fryer, 
Woosley, & Heger 2001). This explosion is called PISN. 

Nozawa et al. (2003) also calculate the dust mass in the ejecta 
of PopIII SNe II and PISNe, carefully treating the radial density 
profile and the temperature evolution. Since it is still debated how 
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Table 1. Summary of grain species. 



Species 


condition" 


Ret* 


density (<5j) 
(g cm- 3) 


C 


u 


1 


2.28 


Si 


u 


2 


2.34 


Si02 


m/u 


3 


2.66 


Fe 


u 


4 


7.95 


FeS 


u 


5 


4.87 


Fe304 


m 


6 


5.25 


AI2O3 


m/u 


7 


4.01 


MgO 


u 


8 


3.59 


MgSiOs 


m/u 


9 


3.20 


Mg2Si04 


m/u 


10 


3.23 



The classifications "m", "u", and "m/u" mean that the species is 
formed in mixed, unmixed, and both supemovae, respectively. 
^ References for optical constants: (1) Edo (1983); (2) Edward (1985); 
(3) Philipp (1985); (4) Lynch & Hunter (1991); (5) Semenov et al. 
(2003); (6) Mukai (1989) (for A < 0.14 ^m, we adopt the values at 
A = 0.14 /im); (7) Toon. Pollack, & Khare (1976); (8) Roessler & 
Huffman (1991); (9) Dorschner et al. (1995); (10) Jager et al. (2003). 

efficiently the mixing of atoms within supemovae (SNe) occurs, 
Nozawa et al. (2003) treat two extreme cases for the mixing of el- 
ements: one is the unmixed case in which the original onion-like 
structure of elements is preserved, and the other is the mixed case 
in which the elements are uniformly mixed within the helium core. 
After examining those two cases, they show that the formed dust 
species depend largely on the mixing of seed elements within SNe, 
because the dominant reactions change depending on the ratio of 
available elements (see Section lTTI for more detailed description). 
Nozawa et al. (2003) predict larger dust mass than Todini & Ferrara 
(2001) for SNe II. How much and what kind of grain species form 
in SNe II and PISNe is still a matter of debate, partly depending 
on the degree of mixing within the He-core, and on the model of 
SNe employed in the calculations; Woosley & Weaver (1995) for 
SNe II in Todini & Ferrara (2001), Heger & Woosley (2002) for 
PISNe in Schneider et al. (2004), and Umeda & Nomoto (2002) for 
PopIII SNe in Nozawa et al. (2003). The formed grain species in 
the calculation of Nozawa et al. (2003) are listed in TableQ 

Some observations have detected infrared radiation from ex- 
tragalactic SNe II (e.g. Dwek et al. 1983; Moseley et al. 1989; 
Kozasa, Hasegawa, & Nomoto 1989). This radiation has been inter- 
preted to be originating from dust formed in SNe II. FIR and sub- 
millimetre (sub-mm) observations of Galactic SN remnants also 
have recently put further constraints on the dust mass formed in 
SNe II (Cas A: e.g. Arendt, Dwek, & Moseley 1999; Dunne et 
al. 2003; Hines et al. 2004; Kepler: e.g. Morgan et al. 2003). Al- 
though FIR and sub-mm observations are useful to know the dust 
amount, the emissivity, which reflects the composition, and the dust 
amount are degenerated in the observed FIR and sub-mm luminos- 
ity. Therefore, in order to constrain the model of dust formation in 
SNe II, another independent information on the dust amount and 
composition is necessary. 

Extinction curves are often used to investigate the dust prop- 
erties (e.g. Mathis 1990). Recently, by using a sample of broad ab- 
sorption line (BAL) quasars, Maiolino et al. (2004a) have shown 
that the extinction properties of the low-redshift (2 < 4, where z is 
the redshift) sample is different from those of the high-z {z > 4.9) 
sample. This result is suggestive of a change in the dust production 



mechanism in the course of galaxy evolution. The highest-redshift 
BAL quasar in their sample, SDSS J104845.05-I-463718.3 (here- 
after SDSS J1048-I-4637) at z = 6.2 shows a red spectrum at the 
restframe wavelength A < 1700 A. (In this paper, all the wave- 
lengths are shown in the restframe of observed galaxies.) However, 
at A > 1700 A, there is no indication of reddening. Then they 
suggest the extinction curve to be flat at A > 1700 A and rising 
at A < 1700 A. Maiolino et al. (2004b) find that the extinction 
curve of SDSS J1048-I-4637 is different from that of low-z BAL 
quasars and is in excellent agreement with the SN II dust models 
by Todini & Ferrara (2001). It is interesting to extend their work 
to various dust formation models in Nozawa et al. (2003). The ex- 
tinction curve should be different from SNe II to PISNe and from 
mixed SNe to unmixed SNe. Dust production in such various con- 
ditions is extensively investigated in Nozawa et al. (2003); hence 
we aim at investigating the extinction curve based on their results. 

In addition to Maiolino et al. (2004a), evidence for dust en- 
richment has been obtained at very high 2 (> 5), where the cosmic 
age is less than 1 Gyr, by the recent sub-mm and millimetre ob- 
servations of distant quasars (Bertoldi et al. 2003; Priddey et al. 
2003). At lower 2 (< 5), direct indications of high-2 dust comes 
from the reddening of background quasars (Fall, Pel, & McMahon 
1989; Zuo et al. 1997). The depletion of heavy elements in quasar 
absorption line systems also supports the presence of dust in distant 
systems (Pettini et al. 1994; Molaro, Vladilo, & Centurion 1998; 
Levshakov et al. 2000; Vladilo 2002; Ledoux, Petitjean, & Sri- 
anand 2003). There are several observations of extinction curves up 
to 2 ~ 1 by taking advantage of the gravitational lensing (Falco et 
al. 1999; and Muiioz et al. 2004). Spectropolarimetric observations 
of two radio galaxies at 2 ^ 1.4 reveal the 2200-A dust feature 
in scattered light (Solorzano-Inarrea et al. 2004). However, most of 
the observations of extinction curves are limited to relatively low 2, 
where dust is not only produced by SNe II, but also by evolved late- 
type stars. In the future, observational samples of extinction curves 
could be extended to high-2 primeval galaxies, where dust is pre- 
dominantly produced in SNe II and/or PISNe. This work could be 
applied to such future observations to reveal the size and composi- 
tion of dust originating from SNe II or PISNe. 

We first describe our theoretical treatment to calculate the ex- 
tinction curves of SN II and PISN dust in Section|2| We examine 
our results, and provide some observationally convenient quantities 
in Section |3| We discuss our results from the observational view- 
point in Section|4] and finally give the conclusion of this paper in 
Section|5| 



2 MODEL 

Our aim is to derive the theoretical extinction curves of dust pro- 
duced in SNe II and PISNe. We adopt the dust production model 
by Nozawa et al. (2003), who investigate various progenitor mass 
of SNe II and PISNe with a careful treatment of physical processes 
(internal mixing, temperature evolution, etc.). 



2.1 Dust production in SNe II and PISNe 

Nozawa et al. (2003) investigate the formation of dust grains in the 
ejecta of PopIII SNe (SNe II and PISNe, whose progenitors are 
initially metal-free). The calculation treats some details compared 
with Todini & Ferrara (2001): (i) the radiative transfer equation in- 
cluding the energy deposition of radioactive elements is solved to 
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calculate the time evolution of gas temperature, which strongly af- 
fects the number density and size of newly formed grains; (ii) the 
radial profile of density of various metals is considered; (iii) un- 
mixed and uniformly mixed cases within the helium core are con- 
sidered. 

In the unmixed case, Nozawa et al. (2003) assume that the 
original onion-like structure of elements is preserved. On the other 
hand, in the mixed case, they uniformly mix all the elements in the 
helium core. They also assume the complete formation of CO and 
SiO molecules, neglecting the destruction of those molecules: no 
carbon-bearing grain condenses in the region of C/0 < 1 and no 
Si-bearing grain, except for oxide grains, condenses in the region 
of Si/O < 1. The formation of CO and SiO may be incomplete be- 
cause of the destruction by energetic electron impact within SNe. 
Todini & Ferrara (2001) treat both formation and destruction of 
CO and SiO, finding that both are mostly destroyed. The decrease 
of CO leads to the formation of carbon grains, which could finally 
be oxidised with available oxygen. The destruction of SiO could 
decrease the formation of grains composed of Si02, MgSiOa, and 
Mg2Si04, and increase other oxidised grains and Si grains. Obser- 
vationally, it is still debated if CO and SiO are efficiently destroyed 
or not. A detailed discussion on this issue can be found in Appendix 
B of Nozawa et al. (2003). 

In the unmixed ejecta, a variety of grain species (Si, Fe, 
Mg2Si04, MgSiOs, MgO, AI2O3, Si02, FeS, and C) condense, 
while in the mixed ejecta, only oxide grains (Si02, MgSiOs, 
Mg2Si04, AI2O3, and Fe304) form. The species are summarised 
in TableQ The difference in the formed species between mixed and 
unmixed cases are mainly derived by the formation of molecules. 
The carbon dust is not produced in their mixed case, because the 
carbon and oxygen are mixed and combined to form CO molecules. 
On the contrary, it forms in unmixed SNe, since there is a carbon- 
rich region at a certain radius of SNe. The formation of SiO 
molecules also affects the formed species: in the mixed ejecta only 
oxide silicate grains form, while in the unmixed ejecta non-oxide 
grains can form in oxygen-poor regions. 

The size of the grains on the location of formation site in the 
ejecta spans a range of 3 orders of magnitude, depending on the 
grain species. The size distribution function summed up over all 
the grain species is approximated by a broken power law. This 
size distribution is different from that of the SN II calculation of 
Todini & Ferrara (2001), which has a typical sizes of 300 A for 
amorphous carbon and 10-20 A for oxide grains. The difference is 
mainly comes from the different treatment of the ejecta: Nozawa 
et al. (2003) consider the density and temperature structures within 
the helium core, while Todini & Ferrara (2001) do not. Schneider et 
al. (2004), based on the model of Todini & Ferrara (2001), find for 
PISNe a large range of dust size, depending on the species. Their 
result is similar to the mixed case in Nozawa et al. (2003). 

We adopt the representative progenitor mass of SNe II as 20 
Mq and that of PISNe as 170 Mq. The size distribution of each 
grain species is almost independent of the progenitor mass, if the 
SN type is fixed (SN II or PISN), except for unmixed SNe II (Sec- 
tion l3.I> . Therefore, we concentrate only those two masses in this 
paper. However, the relative mass ratio among species mildly de- 
pends on the progenitor mass, and we comment on it later showing 
the contribution of each species to the extinction curve (Section 
13. U . We also investigate the mixed and unmixed cases. Therefore, 
we treat four cases: (a) mixed SNe II; (b) unmixed SNe II; (c) mixed 
PISNe; (d) unmixed PISNe, as summarised in Table|2| All the for- 
mulation and the results can be seen in Nozawa et al. (2003). In this 
paper, we assume the grains to be uniform and spherical. 



In FigureQ we show the size distribution adopted in this pa- 
per, where the size distribution function fj{a) is defined so that 
fj (a)da is proportional to the number of grains in a radius interval 
[a, a + da] (j indicates the species). The four figures correspond to 
the four cases in Table|2| The normalisation of fj (a) is discussed in 
Section I7!4| and we only apply an arbitrary normalisation to each 
figure. 

2.2 Calculation of Extinction Curves 

In order to calculate extinction curves, the optical constants for the 
grains are necessary. We adopt the references listed in Table Q for 
the optical constants. By using those optical constants, we calculate 
the absorption and scattering cross sections of homogeneous spher- 
ical grains with various sizes based on the Mie theory (Bohren & 
Huffman 1983). The efficiency factor of extinction, which is de- 
fined as the cross section divided by the geometrical cross section, 
is denoted as Qcxt, j(A, a). This efficiency factor is a function of 
the wavelength A and the dust size a (j denotes the grain species). 

The optical depth of grain j as a function of wavelength, Tx,j, 
is calculated by weighting the cross sections according to the size 
distribution as shown in FigureQ 

r\,j = / vra Qcxt,j(A, a)C/j(a)da, (1) 
Jo 

where C is the normalisation constant related to the column den- 
sity of dust. The determination of C is not important for the extinc- 
tion curve because the extinction curve is presented in the form of 
A\/Av (A)^ is the extinction in units of magnitude at wavelength 
A, and the \/-band wavelength is 0.55 t-tm), where the constant C 
is canceled. The determination of C is necessary when we quantify 
the column density (equation|5}- The extinction in units of magni- 
tude is proportional to the optical depth as 

ylA,j = 1.086rA,i, (2) 

where Ax^j is the extinction of species j in units of magnitude 
as a function of A, and the factor 1.086 comes from 2.51ogjg e. 
The total extinction Ax is calculated by summing Ax,j for all the 
concerning species: 

Ax = ^Ax,j. (3) 
j 



3 RESULTS 

3.1 Extinction curves of various SNe 

In Figure |2[i, we show the extinction curves of the four cases in 
TableU (a) mixed SNe II; (b) unmixed SNe II; (c) mixed PISNe; 
(d) unmixed PISNe. The contribution of each species is also shown. 
The extinction curve is normalised to Av- 

The extinction curves of dust produced by the mixed SNe II 
and PISNe are dominated by Si02 (Figures |2ji and c). Actually, 
Si02 is the most abundant dust component in a mixed 20 Mq SN 
II. However, in a 30 Mq SN II, the production of Mg2Si04 is en- 
hanced by 2.7 times relatively to that of Si02, and in this case, the 
contribution of Mg2Si04 to the extinction curve becomes 2.7 times 
larger As a result, the extinction curve is dominated by the steep 
rising curve of Mg2Si04 for A < 0.14 fxm. For the mixed 170 Mq 
PISN, the extinction curve is also dominated by Si02. Nozawa et 
al. (2003) also examine a larger progenitor mass (200 Mq), where 
Si02 is the dominant species in dust mass. Therefore, we expect 
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Table 2. Models of dust production in supemovae. 



Model 


Progenitor mass 


Mixing 


Rv 


(<Td(V)/md> 


(ad(0.3 ^^m)/m^) 










(10* 


crc? g-i) 


a 


20 


mixed 


2.4 


0.98 


2.8 


b 


20 


unmixed 


3.3 


2.2 


4.4 


c 


170 


mixed 


1.4 


0.75 


3.0 


d 


170 


unmixed 


5.0 


2.1 


4.1 


Galactic'^ 






3.1 


3.4 


6.3 



Quantities derived from observational properties of the Galactic extinction properties (Spitzer 1978; Cardelli et al. 1989; Mathis 1990). We 
should note that it is not necessary to explain the Galactic properties with our models, because the origin, composition, and size of dust are 
different. 





0.0001 0.0010 0.0100 0.1000 1.0000 
a [yLim] 



0.0001 0.0010 0.0100 0.1000 1.0000 
a [n-To] 




Figure 1. Size distribution function of each grain species in (a) the mixed ejecta with the progenitor of 20 Mq, (b) the unmixed ejecta with the progenitor of 
20 Mq, (c) the mixed ejecta with the progenitor of 170 Mq, and (d) the unmixed ejecta with the progenitor of 170 Mq. The correspondence between the 
species and lines is shown in the figures. 



that the extinction curve of mixed PISNe always show a curve sim- 
ilar to the one in Figure|2j:. 

The extinction curve of dust produced by the unmixed SNe II 
(Figure|2j3) is dominated by Mg2Si04 and FeS for A < 0.5 /im. At 
the longer wavelength, (A > 0.5 /im). Si dominates the curve. The 
extinction curves of unmixed SNe II are flatter than those of mixed 



SNe II, because a significant amount of large (a > 0.1 /im) Si 
grains efficiently absorb long-wavelength photons. In a more mas- 
sive SN II with the progenitor mass of 30 Mq, the production of 
Mg2Si04 is enhanced by 2.2 times compared with Si, while the 
amount of FeS is almost the same; consequently the contribution 
of Mg2Si04 to the extinction curve is enhanced by a factor of 2.2 




Figure 2. Extinction curves calculated for the dust production models listed in Table l2l (a) the mixed ejecta with the progenitor of 20 Mq; (b) the unmixed 
ejecta with the progenitor of 20 Mq; (c) the mixed ejecta with the progenitor of 170 AIq; and (d) the unmixed ejecta with the progenitor of 170 Mq. The 
contribution of each species is also shown. Each curve is normalised to the V-band (A = 0.55 fim) value of the total extinction curve. The correspondence 
between the species and lines is shown in the figures. 



relative to the Si contribution, and the extinction curve becomes 
steeper. On the other hand, in less massive unmixed SNe, the pro- 
duction of amorphous carbon is enhanced, leading to the carbon- 
dominated flat extinction curve. Therefore, the extinction curve of 
unmixed SNe II is sensitive to the progenitor mass, and the progen- 
itor mass dependence is shown in Section lT2l 

The extinction curve of unmixed PISNe is also flat, because 
of Si contribution. The curve mildly rises toward ultraviolet (UV) 
because of contribution from Mg2Si04. If the progenitor mass is 
200 Mq, the production of Si is enhanced by twice relative to 
Mg2Si04, and the contribution of Si becomes comparable to that of 
Mg2Si04 even in UV. This makes the extinction curve flatter than 
that shown in Figure|2ll. 

In general, the extinction curves of mixed cases are steeper 
than that of unmixed cases. This is because the dust opacity of the 
mixed SNe II is dominated by small (a ~ 0.01 /im) Si02 grains. 
On the contrary, large Si grains produced in unmixed SNe have a 
large cross section up to the near infrared (NIR). 

All the four extinction curves are compared in Figure|5^. We 
also show the extinction curves of the Galaxy (Cardelli, Clayton, 
& Mathis 1989 with Rv = 3.1) and the Small Magellanic Cloud 



(SMC) (Pel 1992). It is not necessary that our theoretical curves 
explain the Galactic and SMC curves, since in those two environ- 
ments, dust grains are also produced by late-type stars with long 
lifetimes. 

We do not enter deeply the infrared regime. The infrared prop- 
erties of dust grains are addressed in Takeuchi et al. (2004, in prepa- 
ration; see also Takeuchi et al. 2003). 

3.2 Progenitor mass dependence of unmixed SNe II 

As mentioned in Section ITTI the progenitor mass dependence is 
significant in the case of unmixed SNe II. Therefore, we examine 
various progenitor masses examined by Nozawa et al. (2003); i.e. 
the progenitor masses of 13, 20, 25, and 30 Mq. In Figure|3| we 
present the size distribution for each species (Figures|3^, b, c, and 
d correspond to the progenitor masses of 13, 20, 25, and 30 Mq, 
respectively). Based on those size distributions, we calculate the ex- 
tinction curves by the method described in Section lZ2l The extinc- 
tion curve as well as the contribution from each species is shown in 
Figure|4| We see that the 13 Mq extinction curve is dominated by 
carbon grains, whose typical size is large (~ 0.1 /im). Such large 
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Figure 3. Size distribution function of each grain species in the unmixed ejecta with the progenitor masses of (a) 13 Mq, (b) 20 Mq, (c) 25 Mq, and (d) 30 
Mq. The correspondence between the species and lines is shown in the figures. 



carbon grains produce a flat extinction curve as presented in Figure 
|4[i. As the progenitor mass becomes larger, the contribution from 
Si becomes larger. Although the Si extinction curve is flat, other 
species such as Mg2Si04, FeS, Fe, and Si02 contribute to the ris- 
ing curve at short wavelengths in a complex way. As a result, if the 
dust production occurs in the unmixed SNe II, the extinction curve 
tends to be steeper as the progenitor mass becomes larger. The four 
extinction curves are compared in Figurejsjj. 



a factor of ^ 2. Then the astronomical silicate predicts an extinc- 
tion curves with a shallower slope in NIR. However, the difference 
does not affect our results, since the contribution of Mg2Si04 is 
important in UV. The Mg2Si04 data presented in Scott & Duley 
(1996) has a similar optical constants to Draine & Lee (1984) in 
UV and to our adopted values in optical and NIR. 



3.3 Comparison with standard silicate and graphite 

It is useful to compare our prediction with a standard "astronomi- 
cal silicate" in Draine & Lee (1984). Their optical constants in UV 
are based on olivine (Mg, Fe)2Si04 (Huffman & Stapp 1973), and 
the nearest species in our model is Mg2Si04. Our optical constant 
assumed for Mg2Si04 is almost the same as that of "astronomical 
silicate" in UV, but the astronomical silicates have a larger cross 
section in the optical and NIR than our Mg2Si04. We calculate 
the extinction curve of the astronomical silicate by using the op- 
tical constant of Draine & Lee (1984) and the size distribution of 
Mg2Si04 calculated by Nozawa et al. (2003). The extinction curve 
calculated by this method is the same as our extinction curve of 
Mg2Si04 in UV, but at the V band, the difference becomes at most 



The contribution of carbon grains is small in the four cases 
presented in Figure|2| If we use the optical constant of graphite in 
Draine & Lee (1984) and the size distribution of carbon grains in 
Nozawa et al. (2003), we see a weak bump around 2200 A and the 
overall contribution of carbon is reduced by a factor of two. The 
change of carbon optical properties does not affect the total extinc- 
tion. However, Nozawa et al. (2003) show that carbon becomes a 
principal species if the progenitor mass is around 13 Mq. In partic- 
ular, the production of Mg2Si04 is much reduced in a 13 Mq SN 
II. Indeed, carbon grains dominate the extinction curve in 13 Mq 
unmixed SNe. Regardless of which optical constant we adopt for 
carbon grains, we obtain a flat UV extinction curve for a 13 Mq 
unmixed SN II (see Section l4m . 




Figure 4. Extinction curves calculated for the dust production models of unmixed SN II ejecta with the progenitor masses of (a) 13 Mq, (b) 20 Mq, (c) 25 
A/0, and (d) 30 Mq. Each extinction curve is normalised to the value at the V band (A = 0.55 fim) of the total extinction curve. The contribution of each 
species is also shown. The correspondence between the species and lines is shown in the figures. 



3.4 Useful quantities 

The extinction curve is often characterised by the parameter Rv 
defined as 



Ri. 



E{B-V) 



(4) 



where E{B — V) = Ab — Ay, and B indicates the B-band wave- 
length (0.44 nm). Rv roughly quantify the inclination of the ex- 
tinction curve in the optical. The value calculated for each theoret- 
ical extinction curve is shown in Table|2| 

It is also useful to relate the extinction to the dust amount, be- 
cause dust mass also constrain the dust production model in galax- 
ies (e.g. Hirashita & Ferrara 2002). For this aim, we should deter- 
mine the normalisation constant C in equation Q. The constant C 
is determined to realise the total dust column density as 



4 3 

-Tia Sjfj{a)da, 



(5) 



where /i is the gas mass per hydrogen nucleus (assumed to be 1 .4 in 
this paper; Spitzer 1978), mu is the mass of a hydrogen atom, A'^h 
is the column density of hydrogen nuclei, T> is the dust-to-gas mass 



ratio, and 5j is the material density of grain species j. We calculate 
5j from the radius per unit molecule listed in Robie & Waldbaum 
(1968) and the mass of atoms averaged with the isotope ratio in SN 
ejecta. The material density Sj for each species is listed in TableQ 
To relate the dust cross section to the dust mass, we define the 
cross section per unit dust mass, ((Td(A)/md), as 



E 



7ra Qcxt,j(A, a)fj{a)da 



((7d(A)/r7id> = 



E 



4 3 

-Tj-a Sjfj{a)da 



(6) 



Then, the following expression of extinction is derived by using 
equations (T}, 0, 0, and |6}: 



Ax = 1.086AtmHiVHO(crd(A)/md> 



(7) 



We provide the values of {(Jd{V) /nid) (the cross section in the V 
band per unit dust mass) in Table |2| By using this value the dust 
amount can be quantified if we know A^h from other observations 
(for example, observations of Lya absorption line). Then, Ay is 
obtained observationally for primeval galaxies (a system in which 
dust is predominantly supplied by SNe II or PISNe), we obtain the 
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Figure S. Comparison between the extinction curves calculated by our models, (a) Extinction curves of various types of progenitors. The thick solid, dotted, 
dashed, and dot-dashed lines correspond to Models a. b, c, and d in Tablel2l respectively, (b) Extinction curves for various progenitor masses of unmixed Type 
II SNe. The thick solid, dotted, dashed, and dot-dashed lines correspond to the progenitor masses of 13, 20, 25, and 30 Mq, respectively. We also show the 
extinction curves of the Galaxy (dot-dot-dot-dashed line) and the Small Magellanic Cloud (SMC) (thin dashed line) only for the reference. It is not necessary 
that our model explains the Galactic or SMC curve. 



dust-to-gas ratio. The colour excess E{B — V) may be more easily 
obtained, and in this case, Rv hsted in Table |2| could be used to 
derive the extinction Av- 

In the Galactic ISM, iVnMv = 1.9 x 10^^ cm~^ mag"^ and 
2? = 6 X 10"'' (Spitzer 1978). By using equation we obtain 
{od{V) /nid) = 3.4 X 10* cm^ g"^ for the Galactic dust. We have 
obtained similar values in Models b and d (both assume unmixed 
SNe) and significantly smaller ones in Models a and c (both assume 
mixed SNe). 

The dust cross section per mass is also calculated for a UV 
wavelength (0.3 /im). In Table|2| we list ((Td(0.3 /im)/md). The 
difference between the models is relatively small. Therefore, the 
UV dust extinction is a better tracer of the dust column density 
than the optical extinction. 



4 OBSERVATIONAL DISCUSSION 
4.1 Comparison with high-z data 

At present, there are few observational works of the extinction 
curves at 2 > 5, where the cosmic age is shorter than 1 Gyr. 
With this short timescale, the dust is predominantly formed by SNe 
II and PISNe, since their progenitors have short lifetimes while 
evolved low mass stars require longer timescales to evolve and pro- 
duce dust. Maiolino et al. (2004b) use the extinction curve of SDSS 
J1048-I-4637 to test the hypothesis that dust is predominantly sup- 
plied by SNe II. They explain the extinction curve by the dust pro- 
duction model of Todini & Ferrara (2001). They also investigate 
various initial stellar metallicities from to solar after averaging 
the grain properties for different SNe II over the Salpeter stellar 
initial mass function (IMF), and find that their theoretical extinc- 
tion curves agree with the observational data of SDSS J1048-I-4637 
in the whole metallicity range. 

We also use the restframe UV extinction curve of SDSS 
J 1048-1-4637 shown in Figure 2 of Maiolino et al. (2004b). The ex- 
tinction curve is flat at A > 0.17 ^m, and it increases with a smaller 
rate than the SMC extinction curve toward the shorter wavelength 



at A < 0.17 fim. The plausible range derived by Maiolino et al. 
(2004b) is shown by the shaded areas in Figure|6| In Figure|6^, the 
four theoretical curves calculated by our model are shown. The ex- 
tinction curves calculated with the mixed SN models are too steep 
to explain the observational data. 

The models with the unmixed SNe agrees quite well with the 
observational data. The extinction curve of unmixed SNe II de- 
pends on the progenitor mass as shown in Section lT2l In Figurel^J), 
we show the UV extinction curves of unmixed SNe II with the pro- 
genitor masses of 13, 20, 25, and 30 Mq (solid, dotted, dashed, and 
dot-dashed lines, respectively). The flat behaviour of 13 Mq ex- 
tinction curve comes from the carbon contribution, which produces 
a slight bump around 1/A ~ 4 fim~^. The rise toward the shorter 
wavelength is mainly caused by Mg2Si04. Since the production 
of Mg2Si04 is enhanced in the massive progenitors, the extinction 
curve becomes steeper for more massive progenitors. The unmixed 
SN II models are roughly consistent with the current observational 
data at high z. It is interesting to note that the observational ranges 
lie between the flat curve predicted for low-mass (13 A/0 and 20 
Mq) progenitors and the steep curve predicted by high-mass (25 
Mq and 30 A/0) progenitors. Therefore, the mixture of unmixed 
SNe II with various progenitor mass may explain the observational 
extinction curve. 

Therefore, we calculate extinction curves weighted by IMF: 

Axi<j})= / Ax(m)(l,{m)dm, (8) 

where m is the progenitor mass, m\ and are, respectively, the 
lower and upper mass limits of stars which cause SNe II, A\ (m) 
is the extinction calculated for the progenitor mass m (correctly 
weighted for produced dust mass; i.e., A\(m) is large if the pro- 
genitor produces a large amount of dust), and 4>{jn) is the IMF (the 
number of stars with the mass range of [m, m+dm] is proportional 
to (j>{m) dm). The normalisation of 4'{m) is not important in this 
paper, because the extinction curve is always shown after being nor- 
malised by Av and the normalising constant of 4>{m) is cancelled 
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Figure 6. Ultraviolet extinction curves normalised to the extinction at 0.3 /xm. The range observationally derived by Maiolino et al. (2004b) for SDSS 
J1048+4637 is shown by the shaded area in each figure, (a) Model predictions for various supernova models (Models a-d in Table I2I. (b) Results for various 
progenitor masses unmixed SNe 11. The solid, dotted, dashed, and dot-dashed lines correspond to the progenitor masses of 13, 20, 25, and 30 Mq , respectively, 
(c) Averaged extinction curves of unmixed SNe II. The weight of each progenitor mass is determined by initial mass functions, which are pai'ameterised by 
the power-law index ix) and the lower and upper stellai' masses (mi and mu, respectively). The thick solid, dotted, and dashed lines represent the results for 
X = 1.35, 0.35, and 2.35, respectively, with (mi, mu) = (8 Mq, 40 Mq) (Models A, B, and C in Table I3I respectively). The thin solid and dotted lines 
show the results for (mi, mu) = (13 Mq, 40 Mq), and (8 Mq, 30 Mq), respectively, with x = 1.35 (Models D and E in Tablelsl respectively). 



out. We assume the following power-law form of the IMF: 

<j}{m) = i5:m-'"+'> , (9) 

where K is the normalising constant. The Salpeter IMF is repro- 
duced by a; = 1.35. We examine x = 1.35, 0.35, and 2.35. An 
appropriate stellar mass range for SNe II (core collapse SNe) is 
selected as mi = 8 Mq and mu = 40 Mq (Heger & Woosley 
2002). We use the calculated extinction curves of unmixed SNe II 
with m = 13, 20, 25, 30 Mq, and interpolate or extrapolate the 
values to obtain the extinction curves of arbitrary progenitor mass. 
The slopes x — 0.35 (2.35) represents the case where massive 
(less massive) stars are selectively produced. The extinction curves 
weighted for the IMFs are shown in Figure (S];, where the thick 
solid, dotted, and dashed lines represent the results with x — 1.35, 
0.35, and 2.35, respectively. We also investigate the effect of vary- 
ing mi and mu with x = 1.35: the thin solid and dotted lines in Fig- 
ure show the results with (mi, mu) = (13 Mq, 40 Mq) and 



(8 Mq, 30 A^q), respectively. All the examined IMFs are sum- 
marised in Table|2| and are labeled as Models A-E. 

As expected from Figure (SJ), the contribution from massive 
SNe II tends to increase the slope of the extinction curve. Then, 
Models B and D predict steeper extinction curves than the data of 
Maiolino et al. (2004b). Therefore, the contribution from the SNe II 
whose progenitor mass is around 13 Mq is necessary to obtain the 
flat extinction curve consistent with the observational data. How- 
ever, in order to reproduce the rise toward the shorter wavelength 
A < 0.15 fj,m, the contribution from massive SNe II is necessary. In 
particular, Model C predicts an extinction curve slightly inconsis- 
tent around 1/A ~ 8 ^m~^. We should emphasize that the Salpeter 
IMF (Model A) reproduces the observed data very well. 



4.2 Observational strategy 

A large sample of high- 2; quasars and galaxies will be obtained by 
future observations. Since the cosmic age at 2 = 5 is about 1 Gyr, 
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Table 3. Initial mass functions. 



Model 


X 


mi 

(Me) 


mu 
(Mq) 


A 


1.35 


8 


40 


B 


0.35 


8 


40 


C 


2.35 


8 


40 


D 


1.35 


13 


40 


E 


1.35 


8 


30 



the dust supplied by late-type stars does not dominate the total dust 
amount at 2; > 5. Therefore, if future observations collect a large 
spectroscopic sample of quasars at z > 5, we can directly inves- 
tigate the extinction curve of dust, whose source is probably SNe 
II and/or PISNe. Although we have calculated the dust formation 
based on the PopIII progenitors, we can apply our model to metal- 
enriched systems as long as the dust is predominantly formed by 
SNe II and/or PISNe, because the dust composition and size distri- 
bution is much less sensitive to the progenitor metallicity than to 
the progenitor mass. 

The shape of the extinction curve of high-2 galaxies can be 
compared with the theoretical curves in this paper to constrain the 
size and composition of grains. The measurement of extinction en- 
ables us to measure the dust-to-gas ratio of high-z galaxies, if we 
use equation Q and the cross section per dust mass listed in Ta- 
ble|2| Therefore, the dust production history in high-2; universe can 
also be investigated based on this paper. 

Some theoretical works (Bromm & Larson 2004 and refer- 
ences therein) suggest that first stars born from primordial gas 
(PopIII stars) are massive. If the gas metallicity is less than 
IO^'^Zq (Zq is the solar metallicity), massive stars may selectively 
form (Schneider et al. 2003; see also Omukai 2001), causing PISNe 
at the end of their lives. Therefore, the extinction curve of extremely 
metal-poor galaxies could be compared with our curves calculated 
with PISN models. The extinction curve of SDSS J1048-I-4637 
(Maiolino et al. 2004b) has been shown to be fitted by the mod- 
els of SNe n rather than those of PISNe, indicating that SDSS 
J1048-I-4637 is forming stars whose mass is < 30 A/0. This mass 
range is consistent with some metallicity studies of high-z quasars 
(Venkatesan, Schneider, & Ferrara 2003). 

We finally stress that the absorption properties of dust are also 
important for the observation of atoms and molecules, whose de- 
tectability is affected by the dust extinction (e.g. Shibai et al. 2001). 
Therefore, quantifying high-z extinction is crucial to discuss the 
exploration of high- 2: universe by atomic or molecular lines. 



5 CONCLUSION 

We have theoretically investigated the extinction curves of young 
galaxies in which dust is supplied predominantly from Type II su- 
pernovae (SNe II) and/or pair instability supemovae (PISNe). We 
have adopted Nozawa et al. (2003) for compositions and size dis- 
tribution of grains formed in SNe II and PISNe. We have found 
that the extinction curve is quite sensitive to the internal mixing of 
SNe. The extinction curves of mixed SNe II and PISNe are dom- 
inated by Si02 and are charachterised by the steep rise from NIR 
to UV because the main contribution comes from relatively small 
(~ 0.01 fim) grains. On the contrary, the extinction curves of dust 
produced in unmixed SNe II and PISNe are much flatter, because 
of a large contribution from large-sized ~ 0.1 fim Si grains. 



We have also derived the dust cross section per unit dust mass. 
This quantity is useful to estimate the dust column density from 
extinction. The UV extinction trace the dust column density better 
than the optical extinction. The extinction also affects the observ- 
ability of other molecular or atomic lines. The result of this paper 
can be used to estimate the extinction effect in high-z galaxies. 

Finally, our results are compared with a high- 2: extinction 
curve observationally derived for SDSS J1048-I-4637 at z — 6.2 
by Maiolino et al. (2004b). The comparison favours SNe II with- 
out internal mixing as sources of dust grains. The combination of 
various progenitor mass ranging from ~ 10 Mq to ~ 30 Mq ex- 
plains well the observed extinction curve. Our theoretical extinction 
curves could be further utilised when a sample of high-2 extinction 
curves is taken by future observations. 
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